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In this article, we review our recent studies of microwave penetration depth, lower critical fields, 
and quasiparticle conductivity in the superconducting state of Fe-arsenide superconductors. High- 
sensitivity microwave surface impedance measurements of the in-plane penetration depth \ah in 
single crystals of electron-doped PrFeAsOi_j, (y ~ 0.1) and hole-doped Bai_2,Ki:Fe2As2 {x ~ 0.55) 
are presented. In clean crystals of Bai_a:Ka:Fe2As2, as well as in PrFeAsOi-y crystals, the pen- 
etration depth shows flat temperature dependence at low temperatures, indicating that the su- 
perconducting gap opens all over the Fermi surface. The temperature dependence of superfluid 
density A^i,(0)/A^{,(T) in both systems is most consistent with the existence of two different gaps. 
In Bai_i;Ka;Fe2As2, we find that the superfiuid density is sensitive to degrees of disorder inherent 
in the crystals, implying unconventional impurity effect. We also determine the lower critical field 
-ffci in PrFeAsOi_y by using an array of micro-Hall probes. The temperature dependence of Hci 
saturates at low temperatures, fully consistent with the superfluid density determined by microwave 
measurements. The anisotropy of Hci has a weak temperature dependence with smaller values than 
the anisotropy of upper critical flelds at low temperatures, which further supports the multi-gap 
superconductivity in Fe-arsenide systems. The quasiparticle conductivity shows an enhancement in 
the superconducting state, which suggests the reduction of quasiparticle scattering rate due to the 
gap formation below Tc- From these results, we discuss the structure of the superconducting gap in 
these Fe-arsenides, in comparison with the high- Tc cuprate superconductors. 

PACS numbers: 74.25.Nf, 74.25.0p, 74.20.Rp, 74.25. Fy 



I. INTRODUCTION 

Since the discovery of superconductivity in 
LaFeAs(0i_a;F2:) 1], high transition temperatures 
(Tc) up to 56 K have been reported in the doped 
Fe-based oxypnictides 0, i, i, i, [1, 0, i, II- The new 
high-temperature superconductivity in Fe-pnictides 
has attracted tremendous interests both experimen- 
tally and theoretically. The 'mother' materials have 
antifcrromagnetic spin-density-wave order [To| and the 
superconductivity appears by doping charge carriers, 
either electrons or holes. Such carrier doping induced 
superconductivity resembles high- Tc cuprates, but one 
of the most significant differences is the multiband 
electronic structure having electron and hole pockets 
in the Fe-based superconductors. The nature of su- 
perconductivity and the pairing mechanism in such 
systems are fundamental physical problem of crucial 
importance. The first experimental task to this problem 
is to elucidate the superconducting pairing symmetry, 
which is intimately related to the pairing interaction. 

Unconventional superconducting pairings, most no- 
tably the sign-reversin g s+ state, have been suggested 
by several theories 0, Iliriil Q [H, [3 

featuring the 

importance of the nesting between the hole and electron 
bands. This is also in sharp contrast to other multiband 
superconductors such as MgB2, where the coupling be- 
tween the different bands is very weak. Thus the most 
crucial is to clarify the novel multiband nature of super- 
conductivity in this new class of materials. 

In this context, identifying the detailed structure of 



superconducting order parameter, particularly the pres- 
ence or absence of nodes in the gap function, is of pri- 
mary importance. In the electron-doped LnFeAs(0,F) 
or '1111' systems (where Ln is Lanthanoide ions), while 
several experiments 17, Ig] suggest nodes in the gap, the 
tunnelling measurements [19| and angle resolved photoe- 
mission (ARPES) [20j suggest fully gapped superconduc- 
tivity. In the hole-doped Bai„a;Ka;Fe2As2 or '122' sys- 
tem j2l| . ex per imental situation is controversial as well: 
ARPES (2^. |23|. [23 | and lower critical field measurements 
25] support nodeless gaps while /iSR measurements [2^ 
imply the presence of line nodes in the gap function. 

One of the effective ways to judge the presence or ab- 
sence of nodes in the gap is to investigate the proper- 
ties of thermally excited quasiparticles at low tempera- 
tures. Measurements of the London penetration depth 
A are most suited for this, since the quasiparticle den- 
sity is directly related to A(T). In d-wave superconduc- 
tors with line nodes, for example, the low-temperature 
change in the penetration depth 5\{T) — A(T) — A(0) 
shows a linear temperature dependence, as observed in 
high- Tc cuprate superconductors [13]. In contrast, fully 
gapped superconductors without nodes exhibit an expo- 
nential temperature dependence of 5X{T), reflecting the 
thermally activated behavior of quasiparticles. Also no- 
table is that the penetration depth is typically hundreds 
of nm, much longer than lattice parameters, which en- 
ables us to discuss the bulk properties of Fe-based super- 
conductors. This contrasts with other surface sensitive 
techniques such as ARPES and scanning tunneling spec- 
troscopy. 

Another important means to clarify not only the super- 



2 



conducting gap symmetry, but also the multiband nature 
of superconductivity is an accurate determination of the 
lower critical field Hd- In particular, the two-gap su- 
perconductivity in MgB2 manifests itself in the unusual 
temperature dependence of the anisotropy of Hd in the 
superconducting state [H,!!!, IsO]. However, the rehable 
measurement of the lower critical field is a difficult task, 
in particular when strong vortex pinning is present as in 
the case of Fe-arsenides. We also point out that to date 
the reported values of anisotropy parameter strongly vary 
[al spanning from 1.2 (Ref. [sl) up to - 20 



(ReflHlf in 1111 systems. Although the anisotropy pa- 
rameter may be different for different compounds, this 
apparent large spread may be partly due to the effects of 
strong vortex pinning, which lead to large ambiguity in 
some experiments. 

Here we review our microwave surface impedance 
Zs measurements in single crystals of electron- 



doped PrFeAsOi_y 
Bai-xKa;Fe2As2 (x 



(y ^ OA) 
a 0.55) [3|. 



35| and hole-doped 



By using a sensitive 
superconducting cavity resonator, we can measure both 
real and imaginary parts of Zg pricesely in tiny single 
crystals, from which we can extract the in-plane pene- 
tration depth Xab{T) as well as quasiparticle conductiv- 
ity (Ti (T) that yields information on the quasiparticle dy- 
namics. We also measure Hci{T) in PrFeAsOi_j, crystals 
by using an unambiguous method to avoid the difficulty 
associated with pinning (37j . We directly determine the 
field Hp at which first flux penetration occurs from the 
edge of the crystal by measuring the magnetic induction 
just inside and outside the edge of the single crystals, 
with the use of a miniature Hall-sensor array. This al- 
lows us to extract the temperature dependent values of 
the lower critical fields parallel to the c-axis [H^) and to 
the a6-plane {H^^), as well as the anisotropy parameter 
Hd/Hd in single crystals of Fe-based superconductors. 

The penetration depth in Bai_i:Kj;Fe2As2 is found to 
be sensitive to disorder inherent in the crystals. The 
source of disorder may be microscopic inhomogeneous 
content of K, which is reactive with moisture or oxygen. 
The degree of disorder can be quantified by the quasi- 
particle scattering rate 1/r. We find that Xab{T) in crys- 
tals with large scattering rate shows strong temperature 
dependence which mimics that of superconductors with 
nodes in the gap. In the best crystal with the smallest 
l/r, however, \ab{T) shows clear flattening at low tem- 
peratures, giving strong evidence for the nodeless super- 
conductivity. \ab{T) in PrFeAsOi_j^ also shows fiat tem- 
perature dependence at low temperatures. The super- 
fiuid density can be consistently explained by the pres- 
ence of two fully-opened gaps in both hole- and electron- 
doped systems. The lower critical field measurements 
consistently give the saturation behavior of superfluid 
density at low temperatures. Wc find that the anisotropy 
of the penetration depths 7a = Xc/^ab — H^/Hd, where 
Ac and Xab are out-of-plane and in-plane penetration 
depths, respectively, is smaller than the anisotropy of the 
coherence lengths 7^ = £,ab/£.c = H^2/Hc2^ where S^ab and 



are in- and out-of-plane coherence lengths, and 
and are the upper critical fields parallel and perpen- 
dicular to the a6-plane, respectively. This result provides 
further evidence for the multiband nature of the super- 
conductivity. Finally, the quasiparticle conductivity ex- 
hibits a large enhancement in the superconducting state, 
which bears a similarity with high-Tc cuprates and heavy 
fermion superconductors with strong electron scattering 
in the normal state above Tn. 



II. EXPERIMENTAL 
A. Single crystals 

High-quality PrFeAsOi_j^ single crystals were grown 
at National Institute of Advanced Industrial Science and 
Technology (AIST) in Tsukuba by a high-pressure syn- 
thesis method using a belt-type anvil apparatus (Riken 
CAP-07). Powders of PrAs, Fe, Fe2 03 were used as 
the starting materials. PrAs was obtained by react- 
ing Pr chips and As pieces at 500°C for 10 hours, fol- 
lowed by a treatment at 850° C for 5 hours in an evac- 
uated quartz tube. The starting materials were mixed 
at nominal compositions of PrFeAsOo.e and ground in 
an agate mortar in a glove box filled with dry nitrogen 
gas. The mixed powders were pressed into pellets. The 
samples were then grown by heating the pellets in BN 
crucibles under a pressure of about 2 GPa at 1300°C for 
2 hours. Platelet-like single crystals of dimensions up 
150 X 150 X 30 jMu^ were mechanically selected from the 
polycrystalline pellets. The single crystalline nature of 
the samples was checked by Laue X-ray diffraction [38| . 
Our crystals, whose Tc (« 35 K) is lower than the op- 
timum Tc « 51 K of PrFeAsOi_y [3^], are in the un- 
derdoped regime (y 0.1) [40| . which is close to the 
spin-density- wave order [4l|. The sample homogeneity 
was checked by magneto-optical (MO) imaging 37[. The 
crystal exhibits a nearly perfect Meissner state ~ 2 K 
below Tc] no weak links are observed, indicating a good 
homogeneity. The in-plane resistivity is measured by the 
standard four-probe method under magnetic fields up to 
10 T. The electrical contacts were attached by using the 
W deposition technique in a Focused-Ion-Beam system. 

Single crystals of Bai_a;Ka;Fe2As2 were grown at Na- 
tional Institute of Materials Science (NIMS) in Tsukuba 
by a self-flux method using high purity starting materials 
of Ba, K and FeAs. These were placed in a BN crucible, 
sealed in a Mo capsule under Ar atmosphere, heated up 
to 1190°C, and then cooled down at a rate of 4°C/hours, 
followed by a quench at 850°C. Energy dispersive X- 
ray (EDX) analysis reveals the doping level x = 0.55(2) 
[3Q|, which is consistent with the c-axis lattice constant 
c = 1.341(3) nm determined by X-ray diffraction [i^ . 
Bulk superconductivity is characterized by the magneti- 
zation measurements using a commercial magnetometer 
[36| . We find that the superconducting transition tem- 
perature varies slightly from sample to sample [see Ta- 



3 



TABLE I: Properties of the Bai_a;Ka;Fe2As2 crystals we stud- 
ied. In this study, the transition temperature is evaluated 
from the extrapolation of superfluid density Ws — » 0. 



sample 


size (^m"^) 


Tc (K) 


l/r(40 K) (s-i) 


#1 


320 X 500 X 100 


26.4(3) 


27(3) X 10^^ 


#2 


300 X 500 X 80 


25.0(4) 


21(2) X 10^^ 


#3 


100 X 180 X 20 


32.7(2) 


7.8(5) X 10^^ 



blc|I]. Since this is likely related to the microscopic inho- 
mogeneity of K content near the surface, which can be en- 
hanced upon exposure to the air, we carefully cleave both 
sides of the surface of crystal #2 and cut into smaller size 
(crystal #3). For #3, the microwave measurements are 
done with minimal air exposure time. 



B. Microwave surface impedance 




r(K) r(K) 



FIG. 1: Temperature dependence of the surface resistance Rs 
and reactance Xs at 28 GHz in a PrFeAsOi_y single crystal 
(a) and in a Bai_3;K3;Fe2As2 crystal. In the normal state, 
Rs = Xa as expected from Eq. ([1]). The low-temperature er- 
rors in Rs are estimated from run-to- run uncertainties in Q of 
the cavity. Inset in (a) shows the microwave resistivity pi{T) 
(red circles) compared with the dc resistivity in a crystal from 
the same batch (black squares) . The dashed line represents a 
dependence. 



In-plane microwave surface impedance = Rs + , 
where Rg {Xg) is the surface resistance (reactance), is 
measured in the Meissner state by using a 28 GHz TEqu- 
mode superconducting Pb cavity with a high quality fac- 
tor Q ^ 10^ [H, H^]. To measure the surface impedance 
of a small single crystal with high precision, the cavity 
resonator is soaked in the superfluid ^He at 1.6 K and its 
temperature is stabilized within ±1 mK. We place a crys- 
tal in the antinode of the microwave magnetic field H^^ 
(II c axis) so that the shielding current 1^ is excited in the 
ab planes. The inverse of quality factor 1/Q and the shift 
in the resonance frequency are proportional to Rg and 
the change in Xg, respectively. In our frequency range 
Lu /2Tr « 28 GHz, the complex conductivity a — ai ~ i(72 
can be extracted from Zg{T) through the relation valid 
for the so-called skin-depth regime: 

zg=Rg+iXg=(^i^y^\ (1) 

\ai ~ 1(72 J 

In our frequency range, the skin depth Sd is much 
shorter than the sample width, ensuring the skin-depth 
regime. In the superconducting state, the surface reac- 
tance is a direct measure of the superfluid density Ug via 
XgiT) = MowAafc(r) and XjiT) = Aio"s(T)eVm*. In 
the normal state, cti = ne^T/m* ^ a2 gives Rs{T) = 
Xg{T) = {nQUj/iaifl'^ from Eq. IJ), where n is the total 
density of carriers with effective mass m*. Below Tc, the 
real part of conductivity cti is determined by the quasi- 
particle dynamics, and in the simple two- fluid model, 
which is known to be useful in cuprate superconductors 
[27I im , oi is related to the quasiparticle scattering time 
r through ui = {n — ns)e^T/m*(l -t- oj^r^). 



C. Lower critical fields 

The local induction near the surface of the platelet 
crystal has been measured by placing the sample on 
top of a miniature Hall-sensor array tailored in a 
GaAs/AlGaAs heterostructure [1^. Each Hall sensor has 
an active area of 3 x 3 /im^ ; the center-to-center distance 
of neighboring sensors is 20 /im. The local induction at 
the edge of the crystal was detected by the miniature Hall 
sensor located at < 10 /im from the edge. The magnetic 
field Ha is applied for H\\ c and H\\ a6-plane by using 
a low- inductance 2.4 T superconducting magnet with a 
negligibly small remanent field. 



III. RESULTS AND DISCUSSION 
A. Surface impedance 

Figure [1] shows typical temperature dependence of the 
surface resistance Rg and Xg . In the normal state where 
ujT is much smaller than unity, the temperature depen- 
dence of microwave resistivity pi — 2i?^//iow is expected 
to follow p(T) [see Eq. Such a behavior is indeed ob- 
served for PrFeAsOi-j, as shown in the inset of Fig.[lja). 
Below about 100 K pi (T) exhibits a dependence and it 
shows a sharp transition at Tc « 35 K. As shown in Fig.[Tl 
the crystals have low residual Rg values in the low tem- 
perature limit. These results indicate high quality of the 
crystals. We note that the transition in microwave pi{T) 
is intrinsically broader than that in dc p(T), since the 
applied 28-GHz microwave (whose energy corresponds to 
1.3 K) excites additional quasiparticles just below Tc- We 
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FIG. 2: Temperature dependence of 5Xab{T) /S\ab{0) at low 
temperatures, for two crystals of PrFeAsOi-y. The dashed 
lines represent T-linear dependence expected in clean d-wave 
superconductors with line nodes. The solid lines are low-T 
fits to Eq. 121). 



also use Aaf,(0) = 280 nm, which is determined from the 
lower critical field measurements using a micro-array of 
Hall probes in the crystals from the same batch, and 
thus the absolute values of Rg and Xs are determined for 
PrFeAsOi_y. 

For Bai_a;K3,Fe2As2, Zg(T) shows steeper temperature 
dependence in the normal state as shown in Fig. [Ijb). 
This strong temperature dependence above allows 
us to determine precisely the offset of A"s(0)/Xs(40 K), 
since Rg and should be identical in the normal state. 
From this, we are able to determine Aaf,(T)/AQf,(0) and 
ns( T)/ns(0) = A^f,(0)/A^b(r) without any assumptions 
[43| . This also gives us estimates of the normal-state scat- 
tering rate l/r = 1 / ^loTiXlbW = 2tj(X,(T)/X,(0))2, 
which quantifies the degrees of disorder for the samples 
we used [see Table 



B. Penetration depth and superfluid density 

1. PrFeAsOi-y 

In Fig. [5] shown are the normalized change in the in- 
plane penetration depth SXabiT) = \ab{T) — Ao6(0) for 
two crystals of PrFeAsOi_j,. The overall temperature de- 
pendence of 6\ab{T) in these crystals is essentially iden- 
tical, which indicates good reproducibility. It is clear 
from the figure that 5\ab{T) has flat dependence at low 
temperatures. First we compare our data with the expec- 
tations in unconventional superconductors with nodes in 
the gap. In clean superconductors with line nodes, ther- 
mally excited quasiparticles near the gap nodes give rise 
to the T-linear temperature dependence of 5\ab{T) at 
low temperatures, as observed in YBa2Cu307_5 crys- 
tals with d-wave symmetry [13]. In the d-wave case, 
5\ab{T) / \ab{^) ~ ^-^ksT is expected [13, where Aq is 
the maximum of the energy gap A(k). This linear tem- 
perature dependence with an estimation 2Ao/fcsT'c « 4 
\Vt\ [dashed lines in Fig. [2] distinctly deviates from our 
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FIG. 3: Temperature dependence of the superfluid density 
-^a6(0)/Aa6(T) in PrFeAsOi_H crystals. The solid lines are 
the best fit results to the two-gap model [Eq. ((3]l], and the 
dashed and dashed-dotted lines are the single gap results for 
Ai and A2, respectively. Tc is defined by the temperature 
at which the superfiuid density becomes zero. Note that the 
experimental Xs is limited by ^oui5ci/2 above Tc, which gives 
apparent finite values of Us above Tc- 



data. When the impurity scattering rate Pimp becomes 
important in superconductors with line nodes, the in- 
duced residual density of states changes the T-linear de- 
pendence into T^ below a crossover temperature T*^^ de- 
termined by Pimp [Hi- This is also clearly different from 
our data in Fig. [2] If by any chance the T^ dependence 
with a very small slope should not be visible by the exper- 
imental errors below ^ 10 K, then we would require enor- 
mously high T*. However, since no large residual density 
of states is inferred from NMR measurements even in 
polycrystalline samples of La-system with lower Tc [l3| , 
such a possibility is highly unlikely. These results lead us 
to conclude that in contradiction to the presence of nodes 
in the gap, the finite superconducting gap opens up all 
over the Fermi surface. We note that recent penetration 
depth measurements using MHz tunnel-diode oscillators 
in SmFeAs0i_2;Fj^ [ilj and NdFeAsOo.gFo.i [llj are con- 
sistent with our conclusion of a full-gap superconducting 
state in 1111 system. 

In fully gapped superconductors, the quasiparticle ex- 
citation is of an activated type, which gives the exponen- 
tial dependence 



'^A,b(T) 

Aa6(0) 



ttA 

2fcBT 



cxp 



A 



(2) 



at T < Tc/2 [50|. Comparisons between this dependence 
and the low-temperature data in Fig. [2] enable us to es- 
timate the minimum energy Amin required for quasipar- 
ticle excitations at T = K; i.e. Amin/ZcsTc = 1.6 ± 0.1 
for PrFeAsOi_y. 

We can also plot the superfiuid density = 
A^{,(0)/A^^(T) as a function of temperature in Fig. [31 
Again, the low-temperature behavior is quite flat, indi- 
cating a full-gap superconducting state. We note that by 
using the gap Amin obtained above alone, we are unable 
to reproduced satisfactory the whole temperature depen- 
dence of Us [see the dashed lines in Fig. [3] , although a 
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FIG. 4: (color online), (a) Temperature dependence of the 
normalized 28-GHz microwave resistivity pi(r)/pi(40 K) in 
Bai_a;Ka;Fe2As2 single crystals, (b) Normalized superfluid 
density A^(,(0)/A^5(r) for 3 samples with different normal- 
state scattering rates (see Table 



better fit may be obtained by using a larger value of 
A = 1.76kBTc, the BCS value. Since Fe oxypnictides 
have the multi-band electronic structure [5l|, , we also 
try to fit the whole temperature dependence with a sim- 
ple two-gap model 53] 



ns{T) = xnsi(T) + (1 - x)ns2{T). 



(3) 



Here the band 1 (2) has the superfluid density n^i (71.^2) 
which is determined by the gap Ai (A2), and x defines 
the relative weight of each band to ng. The tempera- 
ture dependence of superfluid density nsi(T) {i — 1,2) 
for each band is calculated by assuming the BCS tem- 
perature dependence of superconducting gap A^ (T) [i^ . 
This simple model was successfully used for the two- 
gap s-wave superconductor MgB2 with a large gap ratio 
A2/A1 « 2.6 [30]. Here we fix Ai = Amin obtained from 
the low-temperature fit in Fig. [21 and excellent results of 
the fit are obtained for ^ilksTc — 2.0 and x = 0.6 [see 
Fig. [3] . These exercises suggest that the difference in the 
gap value of each band in this s;i^tem is less substantial 
than the case of MgB2 [H, [H, [Sfl] . 



2. Bax-^K^Fe-iAs-^ 

In Fig.mja) we compare the temperature dependence of 
normal-state microwave resistivity pi = 1 j a\ = 2i?^ / ^^^lo 
[see Eq. ([IJ] for 3 samples of Bai_2:K3;Fe2As2. As men- 
tioned in section 2.1, crystal #3 was cleaved from #2. We 
find that the cleavage dramatically improves the sample 
quality, and crystal #3 exhibits the sharpest supercon- 
ducting transition and the lowest normal-state scattering 
rate l/r [see Table |l]. Figure IHJb) demonstrates the nor- 
malized superfiuid density ns{T)/ns{Q) = Kbi'^) I Kbi^) 
for these 3 samples. We find that crystals with large 
scattering rates exhibit strong temperature dependence 
of superfluid density at low temperatures, which mim- 
ics the power-law temperature dependence of ns{T) in 




FIG. 5: SXabiT) / XabiO) at low temperatures for the cleanest 
crystal #3 of Bai-a;Ki,Fe2As2. The solid line is a fit to Eq. ^ 
with A = i.SOfcsTc. The dashed line represents T-linear de- 
pendence expected in clean d-wave superconductors [13 with 
maximum gap Aq — 2A;sTc [l7| . 



d-wave superconductors with nodes. However, the data 
in cleaner samples show clear flattening at low tempera- 
tures. This systematic change indicates that the super- 
fluid density is quite sensitive to disorder in this system 
and disorder promotes quasiparticle excitations signifl- 
cantly. It is tempting to associate the observed effect 
with unconventional superconductivity with sign rever- 
sal such as the s± state ,11.] . where impurity scattering 
may induce in-gap states in clear contrast to the con- 
ventional s-wave superconductivity [sj, [s^ . Indeed, Tc 
determined by ris — > is noticeably reduced for sam- 
ples with large l/r [Table [I], consistent with theoretical 
studies [5^. At present stage, however, the microscopic 
nature of disorder inherent in our crystals is unclear, and 
a more controlled way of varying degrees of disorder is 
needed for further quantitative understanding of the im- 
purity effects in Fe-based superconductors. In any case, 
our results may account for some of the discrepancies in 
the reports of superfluid density in Fe-arsenides ^2^, ^] . 

For the cleanest sample (#3), the low-temperature 
change in the penetration depth 5Xab{T) = \ab(T) — 
Aaf,(0) is depicted in Fig. [51 which obviously contra- 
dicts the T-linear dependence expected in clean d-wave 
superconductors with line nodes. The low-temperature 
data can rather be fltted to the exponential dependence 
Eq. for full-gap superconductors with a gap value 
A = 1.30fcsrc. This provides compelling evidence that 
the intrinsic gap structure in clean samples has no nodes 
in 122 system. Together with the fact that in fill the 
low-temperature SXab{T) also shows exponential behav- 
ior, we surmise that both electron and hole-doped Fe- 
arsenides are intrinsically full-gap superconductors. 

Various theories have been proposed for the pairing 
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FIG. 6: XlbiO)/Xlb{T) for crystal #3 of Bai-:rK;^Fe2As2 
fitted to the two-gap model Eq. ([3| (blue solid line) with 
Ai = l.nkBTc (dashed line) and A2 = 2A0kBTc (dashed- 
dotted line). Green dotted line is the single-gap fit using 
A = l.SOkBTc- Above Tc, the normal-state skin depth con- 
tribution gives a finite tail. 



symmetry in the doped Fe-based oxypnictides [Til, 
M, Hi, \m, M, HE .58, 59 i^. As confirmed by recent 
ARPES measurements [52j , doped Fe-based oxypnictides 
have hole pockets in the Brillouin zone center and elec- 
tron pockets in the zone edges [51]. It has been suggested 
that the nesting vector between these pockets is impor- 
tant, which favors an extended s-wave order parameter 
(or s± state) having opposite signs between the hole and 
electron pockets fll|, [TJ [l^ [l^.ll5l. [l^ . Our penetration 
depth result of full gap is in good correspondence with 
such an s± state with no nodes in both gaps in these two 
bands. 

As shown in Fig. [SI the overall temperature depen- 
dence of Us in crystal #3 cannot be fully reproduced by 
the single gap calculations. Considering the multiband 
electronic structure in this system [22, 23, 24], we fit the 
data again to the two-gap model Eq. We obtain an 
excellent fit with Ai/fcsTc = 1.17, A^/ksTc = 2.40 and 
X = 0.55. This leads us to conclude the nodeless multi- 
gap superconductivity having at least two different gaps 
in this system. It is noteworthy that the obtained gap 
ratio is comparable to the value A2/A1 w 2 found in 
the ARPES studies '23', ^23] for different bands. A large 
value oi X — 0.55 implies that the Fermi surface with the 
smaller gap Ai has a relatively large volume or carrier 
number, which is also in good correspondence with the 
ARPES results. 



C. Lower critical fields 

In Fig. [7] we show the field dependence of the "local 
magnetization", Modgc = McT^^odgc — Ha, at the edge 
of the crystal, for H\\ c, measured after zero field cool- 
ing. After the initial negative slope corresponding to the 
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FIG. 7: Local magnetization loops for H\\ c, measured by the 
miniature Hall sensor located at < 10 /^m from the edge of 
the crystal. Arrows indicate the field sweep directions. 



Meissner state, vortices enter the sample and Afodgc(-f^a) 
shows a large hysteresis. The shape of the magnetization 
loops (almost symmetric about the horizontal axis) in- 
dicates that the hysteresis mainly arises from bulk flirx 
pinning rather than from the (Bean-Livingston) surface 
barrier [6l[. As shown in Fig. [71 the initial slope of 
the magnetization exhibits a nearly perfect linear depen- 
dence, Modgc — —ctHa- Since the Hall sensor is placed 
on the top surface, with a small but non-vanishing dis- 
tance between the sensor and the crystal, the magnetic 
field leaks around the sample edge with the result that 
the slope a is slightly smaller than unity. 

Figure [5I[a) shows typical curves of B^/^ = ^y^{M + 
aHaY^^ at the edge (circles) and at the center (squares) 
of the crystal, plotted as a function of Ha, the external 
field orientation H\\ c and T = 22 K. The aiJa-term is 
obtained by a least squares fit of the low-field magne- 
tization. The first penetration field Hp corresponds to 
the field ffp(edge), above which B^/"^ increases almost 
linearly, is clearly resolved. In Fig. [5]Ja), we show the 
equivalent curve, measured at the center of the crystal. 
At the center, B^/^ also increases linearly, starting from 
a larger field, ffp(center). 

We have measured the positional dependence of Hp 
and observed that it increases with increasing distance 
from the edge. To examine whether ifp(edge), i.e. Hp 
measured at < 10 /Ltm from the edge, truly corresponds to 
the field of first flux penetration at the boundary of the 
crystal, we have determined the local screening current 
density jodgc = /^(^^(-Bcdgo - -Boutsidc)/Aa; at the crystal 
boundary. Here B^dgc is the local magnetic induction 
measured by the sensor just inside the edge, and Boutsido 
is the induction measured by the neighboring sensor just 
outside the edge. For fields less than the first penetra- 
tion field, jedgc — PHa is the Meissner current, which is 
simply proportional to the applied field {(3 is a constant 
determined by geometry). At Hp, the screening current 
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FIG. 8: (a) Typical curves of y/B (left axis) at the edge (cir- 
cles) and at the center (squares) of the crystal and ^Ajedgc 
(right axis) plotted as a function of Ha for H\\ c at T = 22 K, 
in which Ha is increased after ZFC. (b) The temperature de- 
pendence of the flux penetration fields Hp at the edge and the 
center of the crystal. The insets are schematic illustrations 
of the experimental setup for H\\ c and H\\ afc-plane. (c) 
Temperature dependence of the difference between Hp in the 
center and at the edge (left axis), compared with the rema- 
nent magnetization Mrem (right axis). 



starts to deviate from linearity. Figure[5fa) shows the de- 
viation Ajodgo = Jedge " PHa as a function of Ha- As de- 
picted in Fig.[8ja), y/ Ajcdgc again increases linearly with 
Ha above Hp{edge). This indicates that the Hp{edge) is 
very close to the true field of first flux penetration. 
In Fig. [5]^b) , we compare the temperature dependence 



of ffp(edge) and iJp (center). In the whole tempera- 
ture range, 7?p(center) well exceeds Hp{edge). Moreover, 
iJp (center) increases with decreasing T without any ten- 
dency towards saturation. In sharp contrast, iJp(edge) 
saturates at low temperatures. Figure [8fc) shows the 
difference between Hp measured in the center and at the 
edge, AHp — i7p (center) — _ffp(edge). AHp increases 
steeply with decreasing temperature. Also plotted in 
Fig- IHc) is the remanent magnetization Mrem the 
Ha = value of Modgo on the decreasing field branch), 
measured at near the crystal center. This is proportional 
to the critical current density jc arising from flux pinning. 
The temperature dependence of AHp is very similar to 
that of jc, which indicates that iJp (center) is strongly in- 
fluenced by pinning. Hence, the present results demon- 
strate that the lower critical field value determined by lo- 
cal magnetization measurements carried out at positions 
close to the crystal center is affected by vortex pinning 
effects and might be seriously overestimated. 

The absolute value of Hd is evaluated by taking into 
account the demagnetizing effect. For a platelet sample, 
Hci is given by 



Hci = Hp/ tanh y^OMb/c 



(4) 



where a and b are the width and the thickness of the 
crystal, respectively [6l|. In the situation where H\\ c, 
a = 63 /im and b = 18 fim^ while a = 18 /xm and & = 63 
/im for H\\ a5-plane. These values yield H^^ — 3.22Hp 
and H^i = 1.2477p, respectively. In Fig. [DJa), we plot 
Hci as a function of temperature both for H \\ c and 
H II a6-plane. The solid line in Fig. ^a.) indicates the 
temperature dependence of the superfluid density nor- 
malized by the value at T = K, which is obtained from 
the results in Fig.[3]of a sample from the same batch [sEj. 
Hci{T) is well scaled by the superfluid density, which is 
consistent with fully gapped superconductivity; it does 
not show the unusual behavior reported in Ref. [63l . We 
note that recent JJ^i measurements on NdFeAs(0,F) by 
using Hall probes [6j| also show flat temperature depen- 
dence of Hci, consistent with our results. To estimate 
the in-plane penetration depth at low temperatures, we 
use the approximate single-band London formula. 



fJ-oHci 



^0 



In- 



A 



ab 



Cab 



0.5 



(5) 



where $o is the flux quantum. Using In Aafc/Cat + 0.5 ^ 5, 
we obtain Xab 280 nm. This value is in close corre- 
spondence with the ^SR results in slightly underdoped 
LaFeAs(0,F) [H. 

Next, Fig. [Sl^b) shows the anisotropy of the lower crit- 
ical fields, 7a obtained from the results in Fig. [D^a). 
Here, since the penetration lengths are much larger than 
the coherence lengths for both H\\ ab and H\\ c, the 
logarithmic term in Eq.([5]) does not strongly depend 
on the direction of magnetic field. We thus assumed 
Hci/ Hci — ^c/^ab- The anisotropy 7a ~ 2.5 at very 
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FIG. 9: (a) Lower critical fields as a function of tempera- 
ture in PrFeAsOi-j, single crystals (left axis). The solid line 
(right a:xis) presents the superfluid density ^^bi'^) / ^\b{T) in 
Fig. |3] determined by surface impedance measurements on 
crystals from the same batch, (b) Normalized temperature 
dependence of the anisotropics of Hci (7a , closed circles) and 
Hc2 (7e; closed squares) in FrFeAsOi-y single crystals. The 
anisotropy of Hc2 in NdFeAsO0.82F0.i8 (7e, open squares) 
measured by Y. Jia et a/, [g^ is also plotted. The dashed 
line is a guide to the eye. 



low temperature, and increases gradually with temper- 
ature. In Fig. [D^b), the anisotropy of the upper criti- 
cal fields 7{ is also plotted, where 75 is determined by 
the loci of 10%, 50% and 90% of the normal-state resis- 
tivity measured up to 10 T in a crystal from the same 
batch |32|. Since Hc2 increases rapidly and well exceeds 
10 T just below Tc for H\\ ah, plotting 7^ is restricted 
to a narrow temperature interval. In Fig. [9ljb), we also 
plot the i/c2-anisotropy data from Ref. 66 measured in 
NdFeAsO0.82F0.i8- These indicate that the temperature 
dependence of 7>, is markedly different from that of 7^. 

According to the anisotropic Ginzburg-Landau (GL) 
equation in single-band superconductors, 7a should coin- 
cide with 7j over the whole temperature range. There- 
fore, the large difference between these anisotropics pro- 
vides strong evidence for multiband superconductivity in 
the present system. In a multiband superconductor, 7a 



and 7^ at Tc are given as 



7|(T.) = 7^(Tc).Aii^, (6) 

where (• • •) denotes the average over the Fermi surface, Va 
and Vc are the Fermi velocities parallel and perpendicular 
to the a6-plane, respectively W, '68l|. Vl represents the 
gap anisotropy ((fi^) — 1), which is related to the pair 
potential V{v,v') = Von{v)n{v'). At T = K, the 
anisotropy of the penetration depths is 



7^(0) = (7) 

The gap anisotropy does not enter 7a (0), while 7^ at 
T = K is mainly determined by the gap anisotropy of 
the active band responsible for superconductivity. Thus 
the gradual reduction of 7a with decreasing temperature 
can be accounted for by considering that the contribu- 
tion of the gap anisotropy diminished at low tempera- 
tures. This also implies that the superfluid density along 
the c-axis X^{Q)/X'^{T) has steeper temperature depen- 
dence than that in the plane A^j,(0)/A^^(r). A pro- 
nounced discrepancy between 7^ and 7^ provides strong 
evidence for the multiband nature of superconductivity 
in PrFeAsOi_y, with different gap values in different 
bands. We note that similar differences between 75 (T) 
and 7a(T), as well as A2(0)/A^(T) and Xl,{()) / X^T) , 
have been reported in the two-gap superconductor MgB2 
[11,1331 . We also note that ARPES [H, and Andreev re- 
flection [69| have suggested multiband superconductivity 
with two gap values in Fe-based oxypnictides. 

Band structure calculations for LaFcAsOi-^Fa, yield 
an anisotropy of the resistivity of approximately 15 for 
isotropic scattering 51], which corresponds to 7^ ~ 4. 
This value is close to the observed value. The fact that 7^ 
well exceeds 7^ indicates that the active band for super- 
conductivity is more anisotropic than the passive band. 
According to band structure calculations, there are five 
relevant bands in LaFeAsOi_a;Fa;. Among them, one of 
the three hole bands near the F point and the electron 
bands near the M point are two-dimensional and cylin- 
drical. The other two hole bands near the F point have 
more dispersion along the c axis [5l| , although the shape 
of these Fermi surfaces is sensitive to the position of the 
As atom with respect to the Fe plane, which in turn 
depends on the rare earth [70|. Our results implying 
that the active band is more anisotropic is in good cor- 
respondence with the view that the nesting between the 
cylindrical hole and electron Fermi surfaces is essential 
for superconductivity. This is expected to make these 
two-dimensional bands the active ones, with a large gap, 
and the other more three-dimensional bands passive ones 
with smaller gaps. 
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FIG. 10: (a) Temperature dependence of quasiparticle con- 
ductivity ai normalized by its Tc value at 28 GHz for two 
crystals of PrFeAsOi_i,. The solid line is a BCS calculation 
[zl of ai{T)/ai{Tc) with r = 1.2 x 10"^^ s, which is esti- 
mated from p{Tc) = 77 fi^lcm and Xab{0) = 280 nm. (b) 
(Ti(r) at 28 GHz for crystal #3 of Bai_a;Ka;Fe2As2 normal- 
ized at its 35-K value. The solid line is a BCS calculataion 
with t{Tc) = 4.4 X 10"^^ s. 



D. Quasiparticle conductivity 

Next let us discuss the quasiparticle conductivity 
(Ti(T), which is extracted from Zs{T) through Eq. ([T]). 
The results for PrFeAsOi_j, crystals are demonstrated 
in Fig. [TOTa). Although at low temperatures we have ap- 
preciable errors, it is unmistakable that cri(T) shows a 
large enhancement just below Tc. This enhancement is 
considerably lar ger than the coherence peak expected in 
the BCS theory [nl- 

In Fig. Iior b) we show the temperature dependence 
of the quasiparticle conductivity ai{T)/ai{35 K) in the 
cleanest sample #3 of Bai_a;Ka;Fe2As2. It is evident 
again that below Tc, o'i{T) is enhanced from the normal- 
state values. Near Tc, the effects of coherence factors 
and superconducting fluctuations [27l | are known to en- 
hance ai{T). The former effect, known as a coherence 



■ YBa^CujOg ,5 
O Bai.^K^Fe^Asj #3 




FIG. 11: Quasiparticle scattering rate 1/r as a function of the 
quasiparticle density n„ (T) = n — Us (T) with a comparison 
to the results in YBa2Cu3 06.95 at 34.8 GHz ;4a]. The solid 
and dashed line are fits to the linear and cubic dependence, 
respectively. 



peak, is represented by the solid lines in Fig. [TOl We 
note that in the s± pairing state, the coherence peak in 
the NMR relaxation rate can be suppressed by a partial 
cancellation of total susceptibility x(q) owing to the 
sign change between the hole and electron bands [H, [T^I ■ 
For microwave conductivity, the long wave length limit 
q — > is important and the coherence peak can survive 
[73|, which may explain the bump in cri(T) just below 
Tc. 

At lower temperatures, where the coherence and fluc- 
tuation effects should be vanishing, the conductivity 
shows a further enhancement for the 122 crystal. This 
(Ti(T) enhancement can be attributed to the enhanced 
quasiparticle scattering time r below Tc. The compe- 
tition between increasing t and decreasing quasiparticle 
density n„ (T) — n — Ug (T) makes a peak in cti (T) . This 
behavior is ubiquitous among superconductors having 
strong inelastic scattering in the normal state ^7, 44, 74]- 
Following the pioneering work by Bonn et al. [45|, we 
employ the two-fluid analysis to extract the quasiparticle 
scattering rate 1 /t{T) at low temperatures below ~ 25 K. 
In Fig.fTlTb). the extracted 1/t(T) is plotted against the 
normalized quasiparticle density n„ (T) / n and compared 
with the reported results in the d-wave cuprate supercon- 
ductor YBa2Cu306.95 [lEl- It is found that the scattering 
rate scales almost linearly with the quasiparticle density 
in our 122 system, which is distinct from l/r in cuprates 
that varies more rapidly as ^ n^. Such cubic depen- 
dence in cuprates is consistent with the T^ dependence 
of spin-fluctuation inelastic scattering rate expected in d- 
wave superconcuctors, which have T-linear dependence 
of n„ [791 . In s-wave superconductors without nodes, 
n„(T) and 1/t(T) are both expected to follow exponen- 
tial dependence ^ exp(— A/fc^T) at low temperatures 
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ffE^ . which leads to the linear relation between l/r(T) 
and n„(T). So this newly found relation further supports 
the fully gapped superconductivity in this system. 

IV. CONCLUDING REMARKS 

We have measured the microwave surface impedance 
of PrFeAsOi_j, and Bai_2:K2;Fe2As2 single crystals. We 
also developed an unambiguous method to determine 
lower critical fields by utilizing an array of miniature Hall 
sensors. The penetration depth and the lower critical 
field data both provide saturation behavior of superfiuid 
density at low temperatures for PrFeAsOi_jy crystals. 
The fiat dependence of Xab{T) and ns{T) at low tem- 
peratures demonstrates that the finite superconducting 
gap larger than ~ LGfesTc opens up all over the Fermi 
surface. For Bai_:cK2.Fe2As2 single crystals we find that 
the temperature dependence of the superfiuid density is 
sensitive to disorder and in the cleanest sample it shows 
exponential behavior consistent with fully opened two 
gaps. 

The multiband superconductivity is also supported 
from the anisotropy of iJci, which is found to decrease 
slightly with decreasing temperature. The lower values 
of Hci anisotropy than those of Hc2 suggest that the ac- 
tive band for superconductivity is more anisotropic than 
the passive band. 

The microwave conductivity exhibits an enhancement 



larger than the BCS coherence peak, reminiscent of su- 
perconductors with strong electron scattering. The scat- 
tering rate analysis highlights the difference from the d- 
wave cuprates, which also supports the conclusion that 
the intrinsic order parameter in Fe-As superconductors is 
nodeless all over the Fermi surface. The present results 
impose an important constraint on the order parameter 
symmetry, namely the newly discovered Fe-based high- 
Tc superconductors are fully gapped in contrast to the 
high- Tc cuprate superconductors. 
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